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ABSTRACT 

Highly magnetized, rapidly rotating compact objects are widely argued as central energy sources of 
7-ray bursts (GRBs). After the GRB, such a magnetar-like object may directly lose its rotational energy 
through some magnetically-driven processes, which produce an ultrarelativistic wind dominated possibly 
by the energy flux of electron-positron pairs. The interaction of such a wind with an outward-expanding 
fireball leads to a relativistic wind bubble, being regarded as a relativistic version of the well-studied Crab 
Nebula. We here explore the dynamics of this wind bubble and its emission signatures. We find that 
when the injection energy significantly exceeds the initial energy of the fireball, the bulk Lorentz factor 
of the wind bubble decays more slowly than before, and more importantly, the reverse-shock emission 
could dominate the afterglow emission, which yields a bump in afterglow light curves. In addition, high 
polarization of the bump emission would be expected if a toroidal magnetic field in the shocked wind 
dominates over the random component. 

Subject headings: gamma-rays: bursts — relativity — shock waves — stars: winds, outflows 



1. INTRODUCTION 

The recent observation of high linear polarization dur- 
ing the prompt 7-ray emission of GRB 021206 (Coburn 
& Boggs 2003) suggests that GRBs be driven by highly 
magnetized, rapidly rotating compact objects. Two pop- 
ular scenarios for their birth are the merger of a compact 
binary or the collapse of a massive star (for a recent re- 
view see Meszaros 2002). In both scenarios, a rapidly ro- 
tating black hole surrounded by an accretion disk seems 
to be a common remnant (Narayan, Paczyhski & Piran 
1992; Woosley 1993; Meszaros & Rees 1997a; Paczyhski 
1998). However, a millisecond magnetar has also been 
argued as an alternative interesting product (Usov 1992; 
Duncan & Thompson 1992; Kluzniak & Ruderman 1998; 
Dai & Lu 1998; Spruit 1999; Ruderman, Tao & Kluzniak 
2000; Wheeler et al. 2000). To explain the complex tempo- 
ral feature, the burst itself, in some of these energy models, 
is understood to arise from a series of explosive reconnec- 
tion events in a rising, amplified magnetic field because 
of the Parker instability. This in fact dissipates the dif- 
ferentially rotational energy and magnetic energy of the 
newborn magnetar or accretion disk. 

After the GRB, the remaining object is reasonably as- 
sumed to be a millisecond magnetar or a rapidly rotating 
black hole surrounded by an accretion disk. For the lat- 
ter object, the magnetic field in the disk could have been 
amplified initially by differential rotation to a magnetar- 
like strength of ~ 10^^ G, and particularly, within the 
framework of the coUapsar/hypernova model, such a field 
could be kept, due to longevity (with days or longer) of 
the disk maintained by fallback of the ejecta. During the 
afterglow, the object at the center will directly lose its ro- 
tational energy by the magnetic dipole radiation or the 
Blandford-Znajek mechanism. 

An energy outflow driven magnetically includes three 
components: low-frequency electromagnetic waves, a rela- 
tivistic wind, and a toroidal magnetic field associated with 
the wind. The wind energy flux is unlikely to be baryon- 



dominated, because the initial explosion should have left 
a clean passage with very few baryon contamination for a 
subsequent outflow. The interaction of this outflow with 
an outward-expanding flrcball implies a continuous injec- 
tion of the stellar rotational energy into the fireball. Dai & 
Lu (1998, 2000), Zhang & Meszaros (2001) and Chang, Lee 
& Yi (2002) discussed the evolution of a relativistic fireball 
by assuming a pure electromagnetic-wave energy outflow, 
while Rees & Meszaros (1998), Sari & Meszaros (2000), 
Zhang & Meszaros (2002), and Granot, Nakar & Piran 
(2003) took into account a variable and baryon-dominated 
injection. 

However, based on the successful models of the well- 
observed Crab Nebula (Rees & Gunn 1974; Kennel & 
Coroniti 1984; Begelman & Li 1992; Chevaher 2000), a re- 
alistic, continuous outflow during the afterglow is expected 
to be ultra-relativistic and dominated by the energy flux 
of electron-positron pairs. As in the Crab Nebula, even 
if an outflow from the pulsar is Poynting-flux-dominated 
at small radii, the fluctuating component of the magnetic 
fleld in this outflow can be dissipated by magnetic recon- 
nection and used to accelerate the outflow, which is even- 
tually dominated by the energy flux of e''"e~ pairs within 
a larger radius ~ 10^''' cm (Coroniti 1990; Michel 1994; 
Kirk & Skjseraasan 2003). In the case of an afterglow, 
therefore, it is natural to expect that the central object 
still produces an ultra-relativistic e"'"e~-pair wind, whose 
interaction with the fireball leads to a relativistic wind 
bubble. This can be regarded as a relativistic version of 
the Crab Nebula. 

In this paper, we explore the dynamics of such a wind 
bubble and its emission signatures. In §2 we present ex- 
pressions of the luminosity of a relativistic wind from a 
highly magnetized, rapidly rotating object. In §§3 and 
4 we discuss evolution of the wind bubble and temporal 
features of the radiation respectively. In the final section 
we summarize our findings and give a brief discussion on 
implications of the bubble. 
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2. THE LUMINOSITY OF A RELATIVISTIC WIND 

We assume that a burst itself arises from a series of 
explosive rcconncction events. After the GRB, we are left 
with a highly magnetized, rapidly rotating eompaet objeet. 
Let's first assume that it is a milliseeond magnetar with 
period P, surface magnetic field strength Bs, moment of 
inertia I, radius Rm, and angle between the rotation axis 
and magnetic dipole moment 9. Since such a pulsar loses 
its rotational energy through the magnetic dipole torque, 
the luminosity of a resulting relativistic wind is given by 

L^c^Ax 10^'BI,,RI, ,P-J ergs-i, (1) 

where B±^i4 = sin^/lO^^G, i?M,6 = it!M/10^cm, and 
^ms = ms. Because of spin-down, this luminosity will 
evolve with time as 

T r^^^^ + l^ ^-2 / ~ const., if i < Tm,o, /r,^ 
oc (1 + t/TMfi) I oc if t > Tm,o, 

where t is the observer time in units of day, Tm.o = 
0.58i3^^j4/45i?,i^/gPo jjjs days is the "initial" spin-down 
timescale of the magnetar at the onset of the afterglow, 
Pq = -Po.ms X 1 ms is the rotation period at this time, and 
/45 = //10"''''gcm^. Usov (1992) also assumed that the 
early spin-down could be due to gravitational wave radi- 
ation besides magnetic dipole radiation. However, the lu- 
minosity for gravitational wave radiation depends on the 
stellar ellipticity which is poorly known, and so we ne- 
glected the effect of this mechanism on spin-down. 

We now discuss another case in which the central object 
is a rapidly rotating black hole surrounded by an accretion 
disk. If the amplified magnetic field in the disk does not 
evolve significantly with time during fallback of the ejecta, 
the rotational energy of this black hole will be gradually 
extracted by the Blandford-Znajek mechanism, whose lu- 
minosity is approximated by 

= 1.5 X lO-'^ip|H,i5(^BH/3M0)2a2f(a)crgs-i 
~ 3 X 104^P|H,i5(^BH/3MQ)2(a/0.2)4ergs-i, (3) 

where /(a) = 1 - [(1 ^/l^)/2Y/^ ~ aVS for the 
rotation parameter a <C 1, Mbh is the black- hole mass, 
and Bbh,15 is the disk field strength in units of 10^^ G 
(Lee, Wijers & Brown 2000). We note that a typical value 
(~ lO^^ergs"^) of the wind luminosity in equations (1) 
and (3) has been invoked by Rees & Meszaros (2000) to 
explain the observed iron lines from some GRBs within 
the framework of the coUapsar/hypernova model. Con- 
sidering the rotational energy of the black hole Pbh = 
/(a)MBHc2 ~ (aV8)MBHc2 for a < 1 (Lee et al. 2000), 
and assuming E'bh = L^, we find that the Blandford- 
Znajek mechanism yields spin-down, similarly to the mag- 
netar case, if the accreted angular momentum is neglected 
because its rate seems to be below the torque driven by 
the Blandford-Znajek mechanism at time of days after the 
burst for typical parameters in the coUapsar/hypernova 
model. Thus, we obtain a crude evolution law of the above 
luminosity, 

T r^(^^+l^ ._2 / ~ const., if f < Tbh.o, /.n 

OC (1 + i/TBH,0) I ^ ^-2^ if i > T^^^^^ 



where Pbh.o - l.0B^l^^^{Mjiu/iMQ)-\ao/Q.2)-^ days 
is the "initial" spin-down time of the black hole and oq is 
the "initial" rotation parameter. One can easily see that 
equations (2) and (4) are similar, which implies that our 
discussion of a relativistic wind bubble in the remaining 
text for magnetars should be valid for black holes. 

3. THE BUBBLE DYNAMICS IN THE THIN-SHELL 
APPROXIMATION 

Similarly to the Crab Nebula, a rotating magnetar at 
the center of an afterglow generates a highly relativistic 
wind dominated by the energy flux of e"'"e~ pairs, with 
bulk Lorentz factor of 'y^ 10"^ - lO'^. Atoyan (1999) 
argued that the Crab pulsar initially had 7^, ~ 10* to in- 
terpret the measured radio spectrum of the Crab Nebula. 
We adopt 7^ = lO** as a fiducial value in our calculations. 
Because 7^ is much larger than the Lorentz factor of the 
medium swept up by the fireball, this wind passes through 
a shock front and decelerates to match the expansion ve- 
locity of the swept-up medium. Therefore, a relativistic 
wind bubble, as a result of interaction of the wind with 
the medium, should include two shocks: a reverse shock 
that propagates into the cold wind and a forward shock 
that propagates into the ambient medium. Thus, there are 
four regions separated in the bubble by these shocks: (1) 
the unshocked medium, (2) the forward-shocked medium, 
(3) the reverse-shocked wind gas, and (4) the unshocked 
cold wind, where regions 2 and 3 are separated by a con- 
tact discontinuity. For simplicity, we here assume that two 
initially-forming forward shocks during interactions of the 
fireball both with the medium and with the wind have 
eventually merged to one forward shock, and also neglect 
effects of the baryon loading whose mass is much less than 
the swept-up mass, when the observer time far exceeds the 
initial deceleration timescale. 

We denote rii and P/ as the baryon number density and 
pressure of region "i" in its own rest frame respectively, 
and 7i is the Lorentz factor of region "«" nic;asur(;d in the 
local medium's rest frame. We derive the relative Lorentz 
factor of region 3 measured in the rest frame of region 
4 as 734 ~ (1/2) (7^/73 + 73 /7«-) - 7t«/(273) > 1 for 
7iu S> 73 3> 1, implying a relativistic reverse shock. Be- 
cause the electron number density in the rest frame comov- 
ing with the unshocked wind is rii = Lyj/{'iTrr^j^meC^) 
(where r is the radius of region 3 and me is the electron 
mass), according to the jump conditions for a relativistic 
shock (Blandford & McKee 1976), the pressure of region 
3 is calculated by 

D' ^ 2 2 _ /r\ 

Ps = -^lun.m,c ^ (5) 

Neglecting the presence of the reverse shock and the 
radiative energy loss of region 2, and assuming an ambi- 
ent interstellar medium with constant density of ni, the 
properties of the shocked medium in region 2 should sat- 
isfy the Blandford-McKee adiabatic self-similarity solution 
with the similarity variable at any radius r, 

X = (1 + I67I) (1 - ^) ' (6) 

where 72 = 72 (-R) is the Lorentz factor of the shocked 
medium just behind the forward shock whose radius is 



3 



denoted as i?, and ti = (R/c)[l + 1/(167!)] is the time 
measured in the local medium's rest frame. The radius r 
can thus be expressed as function of x by 



r = R[l 



1 



167: 



1 



X 



1 



I67I 



(7) 



which implies r ~ i? as long as x I672 for an ultrarela- 
tivistic forward shock. This justifies the thin-shell approx- 
imation, in which the width of region 2 is insignificant as 
compared to the shock radius R. 

According to Blandford & McKee (1976), the pressure 
and Lorentz factor of the shocked medium at radius r are 
given by 



72(r-) =72X' 



-1/2 



(8) 
(9) 



where nip is the proton mass. Along the contact disconti- 
nuity, 73 = 72(?') and Pg = ^2(^)1 which yield 



73 = 72X 



X 



-17/12 



-1/2 



(10) 
(11) 



where the thin-shell approximation r R has been con- 
sidered. For an ultrarelativistic, adiabatic forward shock, 
Blandford & McKee (1976) found its total energy. 
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(12) 



In deriving the temporal laws of the similarity variable at 
the location of the contact discontinuity and the Lorentz 
factors of regions 2 and 3, we should note one crucial effect 
that the photons that are radiated from regions 2 and 3 
at the same time in the local medium's rest frame will be 
detected at different observer times. This is because the 
Lorentz factor of region 3 is smaller than 72 (-R) by a fac;- 
tor of x~^^^ so that for a same time interval in the local 
medium's rest frame, i?/c, the emission from region 3 will 
reach the observer at time, 



t ■ 



R 



2^' 



473C 



(13) 



and the emission from region 2 will reach the observer at 

time, 

R 

472C 

Using equations (12) and (14), the Lorentz factor of the 
shocked medium (i.e., region 2) just behind the forward 
shock is found to evolve with time as 



17En 



72 



1/8 



(15) 



From equations (10)- (13), we have the similarity variable 
at the location of the contact discontinuity, 



^ \l7Eo 



-12/17 



= 3.1 



^52 



-12/17 



(16) 



and the Lorentz factor of region 3, 



73 



1/8 



lG247rnimpc5i39/i7 

1/8^-39/136 



(17) 



where L^^47 = L^/IO''"'' erg s ^, E52 = i?o/10'''^ ergs, 



and 



rii and t are in units of 1 cm ^ and 1 day respectively. 
Letting % = 1, we define a critical time 



t„ = 4.9E52L 



'w,i7 



days. 



(18) 



At this time, the injection energy to the fireball signifi- 
cantly exceeds its initial energy. For t < ta, the similarity 

variable x > 1. It should be emphasized that the dy- 
namics denoted by equation (17) is simply calculated by 
equating the pressures of the two-sided shocked fluids at 
the contact discontinuity. In this derivation, we have ne- 
glected any work done on region 2 by region 3 because the 
pressure of region 3 is much less than Pl^iR) at t < tcr- 

Once the observer's time exceeds tcr, the similarity vari- 
able X = 1- -^t this stage, the total kinetic energy of re- 
gion 2 is approximated as -E'kin,2 = {I2 ~ l)-^swC^, where 
-^'-^sw = ('^11 /'A)R^nim.p is the swept-up medium mass. En- 
ergy conservation requires that any increase of kinetic en- 
ergy of region 2 should be equal to work done by region 
3, 

rf£;kin,2 = l^P^dVi, (19) 



where dV^ 



^TTR'^dR' = AT:R^{dRh^) is the volume 
change of region 3 in its own rest frame. Since regions 
2 and 3 should keep velocity equality along the contact 
discontinuity (viz., 72 = 73), we rewrite equation (19) as 



<^72 

dR 



A-KR^[Pi 



(7l 



l)n\mp(? 



272Mswc2 



(20) 



Considering equation (5) and the shock radius R ~ ^'y^ct, 
the solution of equation (20) becomes 



72 = 73 



(- 



1/8 



\ 128'KnimpcH'^ J 



(21) 



where the dependence of 72 on t is consistent with the one 
derived for a pure electromagnetic energy injection by Dai 

& Lu (1998). 

We next discuss the dynamics of a relativistic wind bub- 
ble: In the case of tcr > Tm,o (viz., E52 > 0.46J45^o7ms)> 
the wind bubble should evolve based on equations (15) and 

(17); for tcr < Tm,o (viz., E52 < 0.46/45PoTmJ: however, 
the Lorentz factors of the wind bubble decay initially as 
72 oc t-^/^ and 73 oc f-39/i36 ^ ^ (stage I), subse- 
quently as 72 = 73 oc t"^/"' at t e (tci^Tufi) (stage II), 
and finally again as 72 oc t~^^^ ai t > Tm,o (stage III). 
It should be pointed out that this discussion assumes a 
negligible radiative loss of region 3. However, no matter 
whether region 3 is at stage I or II, and once it enters the 
fast cooling regime at some time, its pressure will begin to 
become much smaller than that of region 2 and then the 
Lorentz factor of region 2 will decay as equation (15). 
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4. SYNCHROTRON RADIATION AND LIGHT CURVES 

In this section we discuss light curves of the emission 
from a rclativistic wind bubble, assuming ta < Tm.o- The 
dynamics above determines the bulk Lorentz factor and 
the thermal Lorentz factors of the accelerated electrons of 
each region as function of time. We consider synchrotron 
radiation from each region at different stages. Accord- 
ing to the standard afterglow model (Mcszaros & Rccs 
1997b; Sari, Piran & Narayan 1998), the spectrum con- 
sists of four power-law segments separated by three break 
frequencies: the self-absorption frequency Va, the charac- 
teristic frequency Vm, and the cooling frequency Vc, with 
the peak flux -Fl/.max- To calculate them, we assume; that 
for region 3 the electron and magnetic field energy den- 
sities are fractions eg and es of the total energy density 
behind the reverse shock (where Ce + £b = 1), but for 
region 2, fractions and S^b of the total energy density 
behind the forward shock (where £,e + S.B < !)• One may 
expect that Cg ^ and/or cb ^ £,b, as suggested in some 
recent studies (Coburn & Boggs 2003; Zhang, Kobayashi 
& Mcszaros 2003; Kumar & Panaitcscu 2003). In addition, 
we assume that the spectral index of the electron energy 
distribution is p for region 3 and q for region 2. 

At stage l{t < tci), the break frequencies and peak fltix 
of region 3 are derived as 

ul^, = 6.8 X 1012^2 g2gV2_^^2^^ 

xLl!li^E-,'/'\Yh'/'^iiz, (22) 

,,I - •? n X inl4.-3/2 r-27/34 

xEU'\^h-''/''Rz, (23) 

pi _ 88^1/2 ^,-1 1-45/34 p-5/68 

-f^i/,max,3 — °°^B,-lIw,4 wA7 ^52 

xn}/^i?Z,^,P/««mJy, (24) 

where = (p - 2)/{p - 1), cb-i = es/O.l, 7^,4 = 
7tu/10*, and £'l,28 is the luminosity distance to the 
source in units of 10^* cm. We here considered 
only the characteristic frequency and the cooling fre- 
quency for the optical to X-ray emission discussed in 
this paper. In our derivation, we used equations (8), 
(10) and (17) to obtain the energy density of region 
3, eg = 4nimpc27|x~-^'^/-^2 = 4nimpc27|x~^^"^^ = 
0.11L^/^^S52'''/'''n?/*t-i9/68ergcm-3, and the mag- 
netic field strength, B'^ = {SireBe'^y/'^ = 
0.52e^/iiLt/l7i?52''/'''n?/'i-i9/i3^ G. We define a cool- 
ing time s through 1/^ 3 = z/J 3 as 

4,3 = 14(5peeeB,-i7-,4)-'*/*^ 

xLzT'ETW^''<i-J^- (25) 

Prom equations (22) and (23), we find that region 3 is in 
the slow-cooling regime for t < 4, 3 but in the fast-cooling 
regime for f > ^qs) in contrast to the standard afterglow 
model (Sari et al. 1998). A larger value of 7^, would imply 
fast cooling in region 3 earlier on, in which case region 2 
could still stay at stage I. For region 2, the break frequen- 
cies and peak fiux are given by 

= 1-6 X 1013^2^2 _i^b'!-i^52'*"'^' Hz, (26) 



ul, = 8.5 X 10'%%lE-,'/'n^H-'/' Hz, (27) 
Fl,m..,2 = 35ej'_,E,2nl/^Dl%, mJy, (28) 

where .g, = {q - 2)/{q - 1), ^e,-i = ^e/0.1, and Cb,-i = 
^b/0.1 (Sari et al. 1998). 

At stage II {tcr <t < Tm,o), since 72 oc we obtain 

the break frequencies and the peak flux for region 3, 

^",3 ^ ^"3 ^ t-\ Fl]^^^^ cx t^'\ (29) 
and for region 2, 

1/^,2 (xt-\ (xt-\ Fl]^^^^^ (X t\ (30) 

We define another cooling time f^s at which 1/^3 = 1/^3 
as follows 

*"3 = 22(5peees,-i7j„,4)"^i'^,47 n^^^^ days. (31) 

Equation (29) is valid only for t < min(tQ^3, Tm.o)- Oth- 
erwise, once region 3 enters the fast-cooling regime or the 
wind luminosity weakens obviously, its pressure becomes 
insignificant as compared to that of region 2 and thus evo- 
lution of the wind bubble comes back to 72 cx fr'^/^ . 

We assume 3 < Tm,o. If tcr < t < tg^g, then region 3 is 
adiabatic and its emission spectrum is determined by equa- 
tion (29). If ^0^3 < t < Tm.o, region 3 becomes fast cooling, 
and the break frequencies and the peak fiux for region 2 
evolve as 2 t'^^'^, v]^2 t^^^'^ and i^^;'^niax,2 A* 
stage III, the emission flux of region 3: ^^^3 oc v~^t~^'^~^^\ 
where /3 is the spectral index at ^0^3 (Kumar & Panaitescu 
2000). 

On the other hand, > ^m.o is assumed. If t^r <t < 
Tm,0) the spectrum and light curves for regions 3 and 2 are 
obtained by using equations (29) and (30), respectively. 
But if t > Tm.o (stage III), because of adiabatic expansion 
of region 3, the break frequencies and the peak flux decay 

as <3 « f-73/48^ ^ ^-73/48 FI%^^, CC ^"47/48 

(Sari & Piran 1999). 

According to the derived scaling laws of the break fre- 
quencies and peak fiux with time at three stages, we can 
obtain the light curve indices for different frequency bands 
in the case of tcr < 2m, o (see Table 1). As an example, 
Figure 1 presents R-band light curves for typical values 
of the model parameters: = 3, Po.ms = 1, 7to,4 = 1, 

Lw,4.7 = 1) P = 9 = 2.5, Ce = 0.9, eB = S.e = S.B = 0.1, 

£■52 = 1, ni = 1 cm~^, and £'l,28 = 1- We can see that the 
emission fiux from region 2 decays rapidly at time < tcr, 
subsequently fades more slowly at time G (fcr, 21vi,o), and 
finally declines based on the initial evolution law (Dai & 
Lu 1998). More importantly, the emission from region 3 
dominates the afterglow emission, which leads to a bump 
in the afterglow light curve. 

5. DISCUSSION AND CONCLUSIONS 

Based on the successful models of the Crab Nebula, we 
discuss the dynamics of a rclativistic wind bubble and its 
emission signatures. Such a wind bubble is naturally ex- 
pected when an ultrarelativistic e+e^-pair wind from a 
highly magnetized, rapidly rotating object at the center of 
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Table 1 

The light curve index a as function of p or g {F„ oc t~"). Definition: Vp = min(i'm, uc) and ;/o = max{i'm, Vc) 



an afterglow interacts with an outward-expanding fireball, 
regardless of whether this object is a millisecond niagnetar 
or a Kerr black hole. We find that when the injection en- 
ergy significantly exceeds the initial energy of the fireball, 
the bulk Lorentz factor of the wind bubble declines more 
slowly than before. In addition, the reverse-shock emission 
could doininate the; afterglow emission, which leads to a 
bump in afterglow light curves. In this paper, we discuss 
the case of t„ < 7m, o- However, even if tcr > 7m, O) a 
bump in Figure 1 still appears for some parameter space 
as we move the late-time light curve of region 3 to early 
times. 

Bump features have been detected in some events (e.g., 
GRBs 970508, 000301C, 021004, and 030329). To in- 
terpret these features, some other models invoked the 
microlensing event (Garnavich, Loeb & Stanek 2000), 
density-jump medium (Dai & Lu 2002; Lazzati et al. 2002; 
Dai & Wu 2003), pure Poynting-flux injection (Dai & 
Lu 1998; Zhang & Meszaros 2001), baryon-dominated in- 
jection (Rees & Meszaros 1998; Sari & Meszaros 2000; 
Zhang & Meszaros 2002; Granot, Nakar & Piran 2003), 
and two-component jet (Berger et al. 2003; Huang et al. 
2003). The magnetic field in the reversely-shocked region 
of the wind bubble seems to consist of two components: 
a toroidal field and a random field. The latter may be 
naturally generated by the relativistic two-stream insta- 
bility (Medvedev & Loeb 1999). If the toroidal field domi- 
nates over the random component, one would expect high 
polarization of the bump emission, which could be used 
to distinguish between our wind-bubble model and other 
explanations. It is well known that the degree of linear 
polarization of synchrotron radiation from the reversely- 
shocked wind is about lisyn = (?'+l)/(p + '''/3) for a large- 
scale toroidal magnetic field. Thus, the polarization of an 



afterglow in our model could be as high as Hgyn ~ 70% 
for p ^ 2 when the reverse shock emission dominates the 
afterglow emission during an obvious bump. Even if the 
total afterglow flux is dominated by the forward shock 
emission, and the reversely-shocked wind provides only a 
small fraction of the total flux, C, then the shocked pair 
wind could still dominate the polarized flux with linear 
polarization of H ~ CHsyn = 7%(C/0.1)(nsy„/0.7). This 
value is larger than the currently observed data of GRBs 
021004 and 030329 (- 2 - 3%, Rol et al. 2003; Lazzati et 
al. 2003; Greiner et al. 2003). If, on the other hand, the 
random- field strength exceeds the toroidal component, the 
calculated degree of linear polarization is significantly less 
than that estimated above. 

We have considered a spherical wind bubble in our pa- 
per. However, there is evidence that GRBs are colli- 
mated into narrow jets, whose kinetic energy is clustered 
at Ej^t ~ 10^° - 10^1 ergs (Frail et al. 2001; Panaitescu & 
Kumar 2002; Berger, Kulkarni & Frail 2003; Bloom, Frafl 
& Kulkarni 2003). On the other hand, relativistic winds 
from millisecond magnetars are roughly isotropic and their 
total energy is 10^^ ergs. Beaming correction gives an 
injection energy of ~ 10^^ ergs within the initial solid angle 
of a jet. This energy is of the same order as Ejef Actually, 
the jet can get more energy from its central magnetar due 
to sideways expansion. This would favor our model. 
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Fig. 1. — R-band light curves of the emissions from regions 2 (dotted line) and 3 (dashed line). The solid line corresponds to the total flux. 



